Effect of nearby Pearl vortices upon the Ic vs B characteristics of planar Josephson 
junctions in thin and narrow superconducting strips 
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In this paper I show how to calculate the effect of a nearby Pearl vortex or antivortex upon the 
critical current Ic{B) when a perpendicular magnetic induction B is applied to a planar Josephson 
junction in a long, thin superconducting strip of width W much less than the Pearl length A = 2A^/d, 
where A is the London penetration depth and d is the thickness (d < A). The theoretical results 
provide a qualitative explanation of unusual features recently observed experimentally by Golod et 
al} in a device with a similar geometry. 

PACS numbers: 74.50. +r,74.78.-w,74.25.-q,74.78.Na 



I. INTRODUCTION 

Golod et al.-*^ recently reported the use of a planar Nb- 
CuNi-Nb Josephson junction of length W = 3.8 to 
detect the presence of a nearby Abrikosov vortex. A hole 
of diameter ~30 nm was fabricated in the Nb film at a 
distance 0.29 /im from the center of the junction. The 
hole could be used to trap a vortex, which carries mag- 
netic flux (/)o = h/2e in the same direction as a positive 
applied magnetic induction B, or an antivortex, which 
carries 4>q in the opposite direction. Without a vortex 
or antivortex in the hole, the Josephson critical current 
Ic{B) vs B exhibited a central maximum with secondary 
peaks roughly resembling the familiar Fraunhofer single- 
slit diffraction pattern. However, when an antivortex was 
trapped in the hole, (i) the central maximum was re- 
placed by a minimum, (ii) the Ic{B) pattern was shifted 
by A<i> K, 00 /2, (iii) an approximate doubling of the pe- 
riodicity appeared on one side of the pattern, leading to 
a clear left-right asymmetry, and (iv) when a vortex was 
trapped in the hole, the Ic{B) pattern was the mirror 
image of that for an antivortex. 

To calculate Ic{B) vs B using the exact geometry and 
material properties used in Ref. 1 would be a very diffi- 
cult numerical problem. Instead, in this paper I consider 
a simpler geometry (see Fig. 1) and solve for Ic{B) vs B 
in the presence of a nearby vortex or antivortex in the 
limit that the strip width W is much less than the Pearl 
length,^ A = 2X^/d, where A is the London penetration 
depth and d is the strip thickness (d < A). This assump- 
tion affords two important simplifications. An applied 
magnetic induction B — zB induces screening currents 
in the film, but when W A, the self-field generated by 
the screening currents can be neglected.^"'* Moreover, in 
this limit a vortex in the strip is best described as a Pearl 
vortex,^ whose properties are totally dominated by the 
1/r sheet-current density circulating around the vortex 
core generated by the gradient of the order parameter's 
phase; within a distance r ^ A from the vortex core the 
vortex's self-field can be neglected. 

Various studies have shown that there is a nonlocal 



FIG. 1: Considered here is a long, thin superconducting strip 
of width W with a planar Josephson junction of width di (bold 
line) at X = separating the two halves a and b. Current leads 
(not shown) symmetrically feed current I to the sample along 
the X direction. A magnetic induction B can be applied in 
the z direction. A possible vortex position on side b is shown 
by the black point with coordinates {xv,yv)- 



relationship between the Josephson-current distribution 
in the vicinity of a Josephson vortex core and the mag- 
netic field these currents generate, and when £ <C A, 
the characteristic length describing the spatial varia- 
tion of the gauge-invariant phase across the junction is 
£ = 4>o/ATrfj.oX^jc, where (/jq = h/2e is the superconduct- 
ing flux quantum and jc (assumed to be independent 
of position) is the maximum Josephson current density 
that can flow as a supercurrent through the junction. 
When ^ A, the characteristic length scale is VIX}^ In 
this paper I assume that the junction length obeys both 
W <^ i and W <C VtK, such that the conditions are 
equivalent to the short-junction limit in sandwich-type 
Josephson junctions. ^^'^^ 

The purpose of this paper is first to review how the 
screening current and the phase gradient induced in re- 
sponse to B affect Ic{B) and then to calculate how Ic{B) 
is affected by the screening current and its phase gradi- 
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ent generated by a vortex or antivortex trapped near the 
junction. 



II. GAUGE-INVARIANT PHASE DIFFERENCE 

In the context of the Ginzburg-Landau (GL) 
theory/^' the superconducting order parameter can be 
expressed as ip = ipofe^'', where tpo is the magnitude of 
the order parameter in equiUbrium, / = |'^|/'^o is the re- 
duced order parameter, and 7 is the phase. The second 
GL equation (in SI units) is 



(1) 



where K = jd is the sheet-current density, A is the vec- 
tor potential, and S = V x A is the magnetic induction. 
Since li" is a gauge- invariant quantity, so is the quantity 
within the parentheses on the right-hand side. Different 
choices for the gauge of the vector potential A result in 
different expressions for 7. 

Consider the planar Josephson junction sketched in 
Fig. 1. With a sinusoidal current-phase relation, the 
Josephson current density in the x direction across the 
junction of width at a; = is Kx{y) = KcSm.A'y{y), 
where Kc = jcd is the maximum Josephson sheet-current 
density and A."/(y) is the gauge-invariant phase difference 
between the left (a) and right (b) superconductors. 



A7(y) =7a(-^,2/) -7b(y,y) 



27r 



di/2 



J-d,/2 



A^{x,y)dx. 

(2) 

I assume here that the induced or applied sheet-current 
densities Xa and Kb on the left- and right-hand sides of 

the jiinction are so weak that the suppression of the mag- 
nitude of the superconducting order parameter is negli- 
gible, such that / = 1. A simple relation between these 
current densities and the gauge-invariant phase difference 
can be obtained by integrating the vector potential A 
around a very narrow rectangular loop of width di in the 
xy plane that just encloses the junction (with the bottom 
end at the origin and the top end at y), neglecting the 
magnetic flux up through the contour, and making use of 
Eq. (1) with / = 1 for those portions of the integration 
along the sides of the junction: 

Aj{y) = A70 + ^ [\Kby{Q,y') - K,y{Q,y')]dy' , 

(3) 

where A70 = ^7(0). In the presence of both an ap- 
plied magnetic induction B and trapped vortices, the 
sheet-current density in general is the vector sum of three 
contributions:^'^^ K = Kj + Kb + K^, where Kj is gen- 
erated by the injection of Josephson currents across the 
junction, Kb is induced by the applied magnetic induc- 
tion B, and K^ is generated by the trapped vortices. 

The short-junction-limit assumption that both W ^ I 
and W <C VIK allows us to neglect the contributions 



from Kj on the right-hand side of Eq. [Z].'^'^^ Thus there 
are only two contribiitions to the sheet-current density 
and gauge- invariant phase difference we need to calculate: 
Kb and A7B induced by the applied magnetic induction 
and Ky and A7^ generated by any nearby trapped vor- 
tices. 



III. A7B INDUCED BY AN APPLIED FIELD 
WHEN NO VORTICES ARE TRAPPED NEARBY 

Let us first calculate the contributions to the sheet- 
current density Kb and the gauge-invariant phase differ- 
ence A7B generated by a perpendicularly applied mag- 
netic induction B = Bz. Since Ksa and the correspond- 
ing phase field ^Ba easily can be obtained by symmetry 
from Ksh and ^Bh, we can calculate only the latter in 
the region a; > and suppress the subscript b. 

With the gauge choice A = —xBy, since V • Kb = 
[sec Eq. (1)], V^jb = must be solved subject to the 
boundary conditions following from KBx{0,y) = and 
KByix,±W/2) = 0, namely 7sx(0,y) = 2'KBy/(j)o and 
^lBy{x,±W/2) = 0, where ^ibx = djB/dx and -/By = 
d^s/dy. The solution for a; > 0, obtained by the method 
of separation of variables, is* (up to a constant) 

, , 8Biy2°2,(_i)"exp[-(2n-M)X]sin[(2n-M)r] 
7B(a;,y) = --2— 2^ 



iBW"' 



71=0 



(2n + 1)3 



e-(^+^^)[-$(-e-2(^+^n, 3,1/2) 
+e2^^$(-e-2(^-*^),3,l/2)], 



(4) 



where X = ttx/W, Y = ny/W, and $(z,s, a) = 
^^gz'^/(A; -I- a)" is the Lerch transcendent.^^ Since 
7b((ii/2,y) in Eq. (2) corresponds to 7B(0,y) and 
la.{-di/2,y) = -7s(0,t/) by symmetry, the gauge- 
invariant phase difference given in Eq. (2) can be ob- 
tained from Eq. (4) as A'ysiy) = — 27^(0, y): 

, , , 16BW^2 °2,(_i)"sin[(2n-M)y] 

Mb {y) = — 2:1— yj — \^ ',.3 

•^^h ^ (2n-M)3 



iBW 

7r2(/)o 

„2iY 



e-*^[$(-e-2*^),3,l/2) 



<I>(-e2«^),3,l/2)]. 



(5) 



The maximum value oiAjsiy) occurs aty = W/2, where 

(6) 



. /TT./^x 14C(3)BTy2 , 
A^BiW/2) = — = 1.705- 



and C(3) = 1.20206 is the Riemann zeta function. Figure 
2(a) shows a plot of A'j'Biy)/AjB{W/2) vs y/iW/2) and 
for comparison sm{TTy /W) vs y/{W/2). If desired, the x 
and y components of the induced sheet-current density 
KB{x,y) can be obtained from Eqs. (1) and (4). 

The Josephson critical current Ic{B), the maximum 
integral of JsTc sin[A7o -|- A7B(y)] over y from —W/2 to 
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FIG. 2: (a) A-ysiy) / A'yB{W/2) from Eq. (5) vs y/{W/2) 
(solid curve). For y < 0, note that AjB{—y) ~ — A7s(y). 
For comparison, the dashed curve shows sin(7r2y/Vl^). (b) 
Ic{B)/Ico from Eq. (5) vs BW'^/<po (solid curve). The dashed 
curve shows the Bessel- function approximation of Eq. (9) . 



W/2, occurs when A70 = ±7r/2, such that 



Ic{B) 



1 

W 



w/2 



W/2 



cos[A7s(y)]rfy 



(7) 



where Ajsiy) is given in Eq. (5) and IcO = KcW . The 
solid curve in Fig. 2(b) shows a plot of Ic{B)/I^ vs 
BW'^/(l)Q. As noted in Ref. 3, the maxima of Ic{B) 
decrease as 1/VB instead of 1/B as in the familiar 
(Praunhofer-like) bulk case. Moreover, the spacings be- 
tween the minima of Ic{B) are not all the same, in con- 
trast to the Fraunhofer pattern. 

Let us define ABi as the value of B at which Ic{B) 
has its first zero, AB2 as the difference of the values at 
which Ic{B) has its second and first zeros, and Ai3„ as 
the difference of the values at which Ic{B) has its nth 
and (n — l)th zeros. For large n, the AB„ approach the 
limiting value 



AS = [7rVl4C(3)](/-o/VK2 ^ l.M2(j)Q/W^ 



(8) 



as pointed out in Refs. 17 and 18. The Ai3„ are smaller 
for small n than for large n. Numerical evaluation of Eq. 
(7) yields the following values for n = 1, 2, 3, 4, and 5: 
ABn/AB = 0.8173, 0.9866, 0.9946, 0.9968, and 0.9979. 
The first minimum of Ic{B) occurs to the left or right of 
the origin i? = at ABi = 1.5O50o/W^^i as can be seen 
in Fig. 2(b). 



If the y dependence of Ajsiv) is approximated by a 
sine function as in the dashed curve in Fig. 2(a), then 
the integral in Eq. (7) can be evaluated in terms of the 
Bessel function Jq with the result 



IciB) 



/14C(3)BI^2 



(9) 



shown as the dashed curve in Fig. 2(b). For large n, the 
spacing between zeros for this approximation to Ic{B) is 
exactly the same as in Eq. (8), but from the well-known 
zeros of Jo{x), we find the following values for n = 1, 2, 
3, 4, and 5: AB^/AB = 0.7655, 0.9916, 0.9975, 0.9988, 
and 0.9993. 



IV. A7„ GENERATED BY A PINNED VORTEX 
WHEN NO MAGNETIC FIELD IS APPLIED 

Let us next calculate the contributions to the sheet- 
current density and the gauge-invariant phase dif- 
ference A7„ generated by a z-oriented Pearl vortex cen- 
tered at {xy,yy) in side b, as shown in Fig. 1. Since 
we are considering the limit W <^ A, we can ignore 
the magnetic field generated by the vortex but we must 
correctly account for the boundary conditions on the 
sheet-current density Ky^ circulating around the vor- 
tex on side b. Because W ^ A, the current density 
on side a is negligibly small {Ky^ = 0). Since we also 
may take the vector potential A to be negligibly small, 
Eq. (1) and V • Ky = yield the equation V^7t., = 0, 
which must be solved subject to the boundary conditions 
following from Ky^{0,y) = and Kyy{x,±W/2) = 0: 
lvx{0, y) = lvy{x, ±W/2) = 0, where 7^^ = d^y/dx and 
Ivy — d-fy/dy. In addition, jy must increase by 27r when 
traversing a closed contour clockwise around the vortex 
axis: V x V7,, = —z2TrS{x — Xy)5{y — yy). The solution, 
obtained using conformal mapping, is 



7v{xv,yv]x,y) = 91n (^^ 



), (10) 



where 9 denotes the imaginary part, ( = x + iy, (y = 
Xy + iyy, and w(C) = ismh{TT(/W). Figure 3 shows 
a plot of the vortex-generated phase when the vortex 
is at {xy,yy) = (0.5W, 0), and Fig. 4 shows a similar 
plot but with more contours for a vortex at (xy,yy) = 
(0.25VF, 0.25VK). If desired, the x and y components of 
the induced sheet-current density KB{x,y) can be ob- 
tained from Eqs. (1) and (10). 

In the limit W <^ A when there is no vortex on side 
a, 7a(-dj/2,i;) = 0, while 7b(rfi/2,y) = -fy{xy,yy;0,y). 
Since we can neglect the vector potential A, Eq. (2) thus 
yields the gauge-invariant phase difference for a Pearl 
vortex (up to a constant) A'yy{y) = — 7^(0;^, y^; 0, y), i.e.. 



A .ir sin(f)-cosh(^)sin(^) ^ 
^""^^^ = -'*^^ i sinh(-)cos(^) J 



(11) 
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FIG. 3: Contour plot of the phase 'yv(xv,yv \ x, y) in the region 
a; > around a Pearl vortex at {xv,yv) ~ (0.5W^, 0). 




FIG. 4: Contour plot of the phase "fv{xv,yv', x, y) in the region 
X > around a Pearl vortex at {xv,yv) = (0.25W, 0.25W). 



The sign of (y) is reversed for a Pearl antivortcx. 

In zero applied field, the Josephson critical current Ic, 
the maximum integral of Kc sinlA-jo + A'yy{y)] over y 
from —W/2 to W/2, occurs in general when tanA7o = 



cos Ajy I sin A7i, , where 



sin A7-U 



cos A7i, 



1 

W 
1 

W 



w/2 

-W/2 
W/2 



iin[A"fy{y)\dy, 



w/2 



cos[Ajy {y)]dy, 



such that, since Ico — KcW, 



Ic = IcoisinA^y + cos A7„ ) 



1/2 



(12) 
(13) 

(14) 



When a singly quantized {N = 1) vortex or antivortcx 
is trapped on the x axis, as in Fig. 3, A7„(y) is an odd 
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FIG. 5: (Color online) Plot of cos[A^v{y)] vs y/W for a singly 
quantized Pearl vortex or antivortex on the x axis (yv/W — 
0) for x„/W = 0.001 (blue), 0.05 (black), 0.175 (red), 0.4 
(purple), and 2 (blue). 
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FIG. 6: Plot of the critical current for a singly (A'' — 1, solid), 
doubly (A^ = 2, dashed), or triply (A^ = 3, dotted) quantized 
Pearl vortex or antivortex at {xv, yv) as a function of its dis- 
tance Xv from the junction for (a) y^ — 0, (b) 0.2W, and (c) 
0.4VF. 



function of y, such that sin A7t, = 0, and 



= |cosA7i,| = 



1-2 tanh — - 
\ W 



(15) 



Thus /c = at one point along the Xv axis, Xy/W = 
tanh"^(l/2)/7r = 0.175. The reason for this behavior is 
illustrated in Fig. 5. Note that |cos A7„| f« 1 for Xy/W = 
0.001 and 2, but that |cosA^| = for Xy/W = 0.175. 
See also Fig. 6(a) for N = 1. 
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Similarly, when a doubly quantized vortex (or antivor- 
tex) is trapped on the x axis and the gauge-invariant 
phase A7„ is doubled, 



Ic_ 

IcO 



1 



4tanh( — 



1 - tanh' 



V W J 



(16) 



such that /c = at two points along the Xv axis, Xy/W = 
0.088 and 0.386 [see Fig. 6(a), = 2]. When a triply 
quantized vortex is trapped on the x axis and the gauge- 
invariant phase Aj^ is tripled. 



Ic_ 
IcO 



1-2 tanh 



V W J 



3-8tanh2f^ 
W 



6 tanh 



TTXy 

w 



(17) 



such that /c = at three points along the 
Xy/W = 0.059, 0.232, and 0.513 [see Fig. 6(a), N = 3]. 
However, when the vortex is trapped at a position off 
the X axis, as in Fig. 4, the zeros of Ic are replaced by 
minima. As shown in Fig. 6(b) and (c), Ic vs Xy exhibits 
one minimum for a singly quantized vortex {N = 1), two 
for a doubly quantized vortex (N — 2), and three for a 
triply quantized vortex (N — 3). 



V. A7 GENERATED WHEN A MAGNETIC 
FIELD IS APPLIED IN THE PRESENCE OF A 
PINNED VORTEX OR ANTIVORTEX 

We are now in a position to calculate how the Ic{B) 
characteristics calculated in Sec. Ill are affected by the 
presence of a vortex, described in Sec. IV. Since the 
resulting gauge-invariant phase difference A7 is (aside 
from a constant) simply the sum of the contributions 
A7S and A7^ (and the resulting sheet-current density K 
is the sum of Kb and Ky)^ the junction critical current 
in the presence of both an applied magnetic induction 
B = zB and a z-directed Pearl vortex at {xy, yy) is given 
by 



Ic (B) //co = (sin A7 +COSA7 



a/2 



(18) 



where Icq = KcW and the averages are calculated as in 
Eqs. (12) and (13) but with A-yiy) = A-ysiv) + A7„(j/). 
For the case of an applied magnetic induction and an 
antivortex, the sign of A'yy{y) is reversed, and the av- 
erages are calculated with Aj{y) = A^B{y) — Ajy{y). 
For a vortex or an antivortex on the x axis, since both 
^7s(y) s-nd A'yy{y) when yy — are odd functions of y, 
sin A7 = 0, and Ic{B)/Ico = |cosA7|. 

Figure 7 exhibits the interesting behavior of how the 
pattern of Ic{B) vs B depends upon the position of a 
nearby Pearl antivortex at x = Xy and y = 0. First let 
us focus on the behavior of /c(0), which corresponds to 
the case discussed in Sec. IV. In the limit /w 
[blue point in Fig. 7(a)] the critical current is simply 
IcO = KcW . However, as the antivortex moves away from 
the junction, /c(0) drops to smaller values [black point at 
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FIG. 7: (Color online) Behavior of the critical current for an 
antivortex at x = a:„ and y — = 0: (a) Normalized critical 
current Ic/IcO vs Xv/W a.t B — 0. (b) Normalized position 
of the minimum of Ic{B) nearest the origin, BminW^ /cj>o vs 
x^/W. (c) Normalized critical current Ic{B)/Ico vs BW^/(f>o 
for Xv/W = (blue dotted curve), 0.05 (black solid curve), 
and 0.175 (red dashed curve, (d) Normalized critical current 
Ic{B)/ho vs SWV^o for x^/W = 0.175 (red dashed curve), 
0.4 (purple solid curve), and 00 (blue dotted curve). 



Xy/W = 0.05 in Fig. 7(a)] and becomes zero [red point 
at Xy/W = 0.175 in Fig. 7(a)]. As Xy increases further, 
/c(0) rises [purple point at Xy/W = 0.4 in Fig. 7(a)] and 
approaches Ico = KcW as Xy/W — >■ 00 as shown in Fig. 
7(a). 

Next let us focus on the minimum at _B = Bmim 
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where initially in the limit Xy/W [blue point at 
Xy/W = in Fig. 7(b) and blue dotted curve in Fig. 
7(c)] BrmnW^/^a = -1.505, to the left of the origin of 
Fig. 7(c). As the antivortex moves away from the junc- 
tion, this minimum moves to the right, as shown by the 
black point at Xy/W = 0.05 in Fig. 7(b) and the black 
solid curve in Fig. 7(c). When the antivortex reaches 
the point Xy/W = 0.175 [red point in Fig. 7(b) and red 
dashed curve in Fig. 7(c) and (d)], the minimum occurs at 
B = 0. As Xy increases further, this minimum continues 
to move to the right [purple point at Xy/W = 0.4 in Fig. 
7(b) and purple solid curve in Fig. 7(d)]. As Xy/W oo, 
the minimum occurs at BminW^ / (po = +1.505, as shown 
in Fig. 7(b), and the pattern of Ic{B)/Ico becomes the 
same as in the absence of a vortex or antivortex [blue dot- 
ted curve in Fig. 7(d) or black solid curve in Fig. 2(b)]. 

Now let us examine the behavior of the pattern of 
primary and secondary maxima in Ic{B)/Ico as Xy in- 
creases. When Xy/W — 0, there is a secondary maxi- 
mum of height 0.34 at BW^ /ipo = —2.33 and a primary 
maximum of height 1.00 at BW'^/(j)o = [blue dotted 
curve in Fig. 7(c)]. As Xy/W increases, the secondary 
maximum movc;s to the right and grows in height until it 
becomes the primary maximum, while the primary maxi- 
mum also moves to the right but decreases in height until 
it becomes a secondary maximum. These changes can be 
seen in Fig. 7 in the progression of the black solid curve 
in (c), the red dashed curve in (c) and (d), the purple 
solid curve in (d), and the blue dotted curve in (d). 

Finally we note the asymmetry of the pattern of the 
primary and secondary maxima as Xy increases. Al- 
though Ic{B)/Ico has mirror symmetry in the limits 
Xy/W — > and Xy/W —J- oo, as shown by the blue dot- 
ted curves in Fig. 7(c) and (d), this symmetry is bro- 
ken for intermediate values of Xy/W. The asymmetry is 
most pronounced for < Xy/W < 0.175. For example, 
for Xy/W = 0.05 [black solid curve in Fig. 7(c)], the 
secondary maxima decrease monotonically for increas- 
ing B > and have nearly the same period. How- 
ever, for negative B the heights of the secondary max- 
ima are irregular, and the larger secondary maxima have 
approximately double the period of those for S > 0. 
These effects also occur for Xy/W = 0.175 [red dashed 
curve in Fig. 7(c) and (d)] but are less pronounced. For 
Xy/W = 0.4 [purple curve in Fig. 7(d)], asymmetry is 
present, with the secondary peaks for B < lower than 
those for B > 0, but the period doubling is no longer 
present. 

The basic reason for the asymmetry of Ic{B) about 
B = when j/^ = is that Ic{B)/Ico = IC0SA7], where 
A7(y) = Ajsiy) - ^Iviu)- Ahhough A7b(?/) is an an- 
tisymmetric function of B, — A7„(y) is independent of 
B, so that when B ^ 0, A7(y) is neither symmetric nor 
antisymmetric about B = 0. The effects of this asym- 
metry can be very pronounced, as seen in the example 
of the black solid curve in Fig. 7(c). Plots of Aj{y) and 
cos[A7(r/)] calculated for Xy/W = 0.05, shown in Fig. 8 
at the first three maxima for i? > and Fig. 9 at the 



< 



10 
5 


-5 

-JO 



(a) 


BW^I<l>o = 4.373. ■ 


~~ 0.332 









-0.4 



-0.2 0.0 

ylW 



0.2 



0.4 



1.0 
0.5 



^ 0.0 



-0.5 
-1.0 



\ k 


\ i \: 2.506-- 


\ ' /'■ 
\ ' / '.' 

\ / ',' 


'A': ! \ = 4.373 


7; 









-0.4 -0.2 0.0 0.2 0.4 
y/W 

FIG. 8: (a) A-yiy) and (b) cos[A7(;y)] when x^/W = 0.05 
and yv — [black solid curve in Fig. 7(c)] for the first three 
maxima for positive B, for which BW^ /cl>o and Ic{B)/Ico are 
0.332 and 0.759; 2.506 and 0.323; and 4.373 and 0.236, re- 
spectively. 
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FIG. 9: (a) A-fiy) and (b) cos[A7(y)] when x^/W = 0.05 
and yv ~ [black solid curve in Fig. 7(c)] for the first three 
maxima for negative B, for which BW^/(j)o and Ic{B)/Ico 
are -2.069 and 0.750; -4.087 and 0.013; and -5.797 and 0.363, 
respectively. 
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FIG. 10: Theoretically calculated normalized critical current 

Ic{B)/Ico vs BW"^ /(t>o in the absence of a vortex or antivor- 
tex (dashed curve) and in the presence of an antivortex at 
Xv/W = 0.077 and j/ = j/„ = (solid curve). 



FIG. 11: Theoretically calculated normalized critical current 

Ic{B)/IcO vs BW"^ /(f)a in the presence of a doubly quantized 
antivortex {N = 2, dashed) and a triply quantized antivortex 
{N = 3, solid) at Xy/W = 0.077 and y = = 0. 



first thrc'c maxima for B < 0, show the dramatic differ- 
ences responsible for the asymmetry of Ic{B) / Ico and the 
approximate period doubling for B < 0. 

Numerical calculations of how the Ic{B) vs B patterns 
for a doubly {N — 2) or triply {N — 2) quantized an- 
tivortex on the X axis evolve as Xy increases from zero 
to values of order W or larger reveal behavior similar to 
those for a singly (TV = 1) quantized antivortex shown in 
Fig. 7. As Xy increases, the patterns shift to the right, 
and for intermediate values of Xy the maxima decrease 
monotonically for B > but have irregular heights for 
B < 0. The chief difference from the behavior for iV = 1 
is that /c(0) passes through zero twice for = 2 and 
three times for iV = 3 in accordance with Fig. 6(a). 

VI. SUMMARY 

In this paper I first reviewed how the gauge-invariant 
phase difference A'y{y) across a planar Joscphson junc- 
tion in a long, thin superconducting film is affected by the 
sheet-current distributions on opposite sides of the junc- 
tion. The assumptions that T4^ <C A and W £ made it 
possible to calculate; the two relevant independent contri- 
butions to the gauge- invariant phase difference: AjB{y) 
due to the perpendicular applied magnetic induction B 
and Ajy [y) due to a nearby trapped Pearl vortex or an- 
tivortex. After calculating the critical current Ic of the 
junction for these two contributions separately, I calcu- 
lated Ic{B) when both B is applied and a vortex or an- 
tivortex is near the junction. 

The features observed in the calculated Ic{B) vs B 
characteristics show many features in common with the 
experimental Ic{B) vs B characteristics observed recently 
by Golod et al.^ The dashed curve in Fig. 10 shows the 
calculated curve of Ic{B) vs B in the absence of a vor- 
tex or antivortex, which shows perfect mirror symmetry 
about i? = 0, a primary maximum at S = and sec- 
ondary maxima of monotonically decreasing heights for 
increasing \B\. The corresponding experimental curve 
[Fig. 3(a) in Ref. 1] shows approximate mirror symmetry 



about i? = 0, a primary maximum at B = and sec- 
ondary maxima, which generally decrease in height for 
increasing \B\ but not monotonically. The spacings of 
the minima and maxima along the B axis increase for 
increasing \B\, as expected from the discussion in the 
paragraph containing Eq. (8). 

Figure 3(b) in Ref. 1 showed Ic{B) vs B for one an- 
tivortex trapped in a hole at a distance 0.29 /zm from the 
center of the junction, whose length was 3.8 /um. In our 
model of the experiment, this corresponds to having an 
antivortex at Xy/W = 0.29/3.8 = 0.077 and yy = 0, for 
which the calculated Ic{B) vs B is shown by the solid 
curve in Fig. 10. This curve shows a primary maximum 
shifted to the left of B = with secondary maxima to the 
right of the primary maximum monotonically decreas- 
ing for increasing B, and secondary maxima of irregu- 
lar heights to the left of the primary maximum showing 
an approximate period doubling. The experimental plot 
of Ic{B) vs B, shown in Fig. 3(b) in Ref. 1, exhibits 
similar features: a primary maximum shifted to the left 
of i? = 0, secondary maxima to the right of the pri- 
mary maximum monotonically decreasing for increasing 
B, and secondary maxima of irregular heights to the left 
of the primary maximum showing an approximate pe- 
riod doubling. However, the experimental Ic{B) showed 
a minimum at i? = 0, while the theoretical curve has this 
minimum shifted to the left of the origin. 

Figure 11 shows theoretical predictions of Ic{B) vs 
B patterns for a doubly quantized antivortex {N = 2, 
dashed) or a triply quantized antivortex {N = 3, solid) 
trapped in the hole at Xy/W = 0.077 and yy = 0. Note 
the monotonic decrease of the maxima for B > and 
irregular heights for B < 0. 

Although the theoretical model docs not assume the 
exact geometry and material properties of the sample 
used in Ref. 1, the theoretical results presented here pro- 
vide a good qualitative and semi-quantitative explana- 
tion of the experimental results. 
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